Correlated perturbations from inflation and the cosmic microwave background 
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We compare the latest cosmic microwave background data with theoretical predictions including 
correlated adiabatic and CDM isocurvature perturbations with a simple power-law dependence. 
We find that there is a degeneracy between the amplitude of correlated isocurvature perturbations 
and the spectral tilt. A negative (red) tilt is found to be compatible with a larger isocurvature 
contribution. Estimates of the baryon and CDM densities are found to be almost independent of 
the isocurvature amplitude. The main result is that current microwave background data do not 
exclude a dominant contribution from CDM isocurvature fluctuations on large scales. 
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Increasingly accurate measurements of temperature 
anisotropies in the cosmic microwave background sky of- 
fer the prospect of precise determinations of both cos- 
mological parameters and the nature of the primordial 
perturbation spectra. The recent Boomerang [y, DASI 
and Maxima [g| data have shown evidence for three 
peaks in the cosmic microwave background (CMB) tem- 
perature anisotropy power spectrum as expected in in- 
flationary scenarios. In this context the CMB data sup- 
port the current 'concordance' model based on a spatially 
flat Friedmann-Robertson- Walker universe dominated by 
cold dark matter and a cosmological constant Q . In ad- 
dition, the CMB data no longer shows any signs of being 
in conflict with the big bang nucleosynthesis data ^ . 

In the studies which have estimated the cosmological 
and primordial parameters with these new data sets, only 
the case of purely adiabatic perturbations has been con- 
sidered so far. That is, the perturbation in the relative 
number densities of different particle species is taken to 
be zero. Although this assumption is justified for pertur- 
bations originating from single field inflationary models, 
it does not necessarily follow when there is more than 
one field present during inflation Iq-lO . Other possible 
primordial modes are isocurvature |1 1|, 12| (also referred 
to as "entropy") modes in which the particle ratios are 
perturbed but the total energy density is unperturbed in 
the comoving gauge. 

Most previous studies have examined the extent to 
which a statistically independent isocurvature contribu- 
tion to the primordial perturbations may be constrained 
by CMB and large-scale structure data ||l^,Q. It has 
recently been shown that multi-fleld inflationary models 
in general produce correlated adiabatic and isocurvature 
perturbations [[7|-pl|. These correlations can dramati- 
cally change the observational effect of adding isocurva- 
ture perturbations llqJl2]. Up until now, only the case 
of scale-invariant correlated adiabatic and entropy per- 
turbations has been considered. Trotta et al. ilql found 



(with an earlier CMB dataset) that in this case the cold 
dark matter (CDM) isocurvature mode was likely to be 
very small if not entirely absent, though they did find 
that a neutrino isocurvature mode contribution |12(| was 
not ruled out. In this letter we examine to what extent 
a correlated CDM isocurvature mode is consistent with 
the recent CMB data when a tilted power law spectrum 
is allowed. 

Non-adiabatic perturbations are produced during a 
period of slow-roll inflation in the presence of two or 
more light scalar fields, whose effective masses are less 
than the Hubble rate. On sub-horizon scales, fluctua- 
tions remain in their vacuum state so that when fluctua- 
tions reach the horizon scale their amplitude is given by 
54>^^ ~ {H^l 2t:) di where the subscript * denotes horizon- 
crossing and di are independent normalised Gaussian 
random variables, obeying (didj) — Sij. The total co- 
moving curvature and entropy perturbation at any time 
during two-field inflation can quite generally be given in 
terms of the field perturbations, along and orthogonal to 
the background trajectory, as B 
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where 9 is the angle of the inflaton trajectory in field 
space. Although the curvature and entropy perturbations 
are uncorrelated at horizon-crossing, any change in the 
angle of the trajectory, 6, will begin to introduce correla- 
tions B. Further correlations may be introduced by the 
model dependent dynamics when inflation ends and the 
flelds' energy is transformed into radiation and/or dark 
matter. The comoving curvature perturbation, T^rad, 
on large-scales during the radiation-dominated era is re- 
lated to the conformal Newtonian metric perturbation, 
<&, by 7?.rad — 3<f>/2. The isocurvature perturbation is 
^rad = <5/0cdm/Pcdm — i^/ ^)Sp-f / P-y ^^^d remains constant 
on large scales until it re-enters the horizon. On large 
scales the CMB temperature perturbation can be ex- 



pressed in terms of the primordial perturbations 
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The general transformation of linear curvature and en- 
tropy perturbations from horizon-crossing during infla- 
tion to the beginning of the radiation era will be of the 
form 
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Two of the matrix coefficients, T-jm = 1 and Tsti = 0, 
are determined by the physical requirement that the 
curvature perturbation is conserved for purely adiabatic 
perturbations and that adiabatic perturbations cannot 
source entropy perturbations on large scales |17| . The re- 
maining terms will be model dependent. If the fields and 
their decay products completely thermalize after inflation 
then Tss = and there can be no entropy perturbation 
if all species are in thermal equilibrium characterised by 
a single temperature, T. This means that it is unlikely 
that a neutrino isocurvature perturbation could be pro- 
duced by inflation unless the reheat temperature is close 
to that at neutrino decoupling shortly before primordial 
nucleosynthesis takes place. On the other hand, a cold 
dark matter species could remain decoupled at temper- 
atures close to, or above, the supersymmetry breaking 
scale yielding T55. The simplest assumption being that 
one of the fields can itself be identified with the cold dark 
matter m. 

The slow evolution (relative to the Hubble rate) of light 
fields after horizon-crossing translates into a weak scale 
dependence of both the initial amplitude of the pertur- 
bations at horizon crossing, and the transfer coefficients 
T-jis and T^s- Parameterising each of these by sim- 
ple power-laws over the scales of interest, requires three 
power-laws to describe the scale-dependence in the most 
general adiabatic and isocurvature perturbations. 
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The generic power-law spectrum of adiabatic perturba- 
tions from single field inflation can be described by two 
parameters, the amplitude and tilt, A and n. Uncor- 
related isocurvature perturbations require a further two 
parameters, whereas we now have in general six param- 
eters. The dimensionless cross-correlation 
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is in general scale-dependent. 

We will investigate in this letter the restricted case 
where all the spectra share the same spectral index and 
hence A is scale-independent. This might naturally arise 



in the case of almost massless fields where the scale- 
dependence of the field perturbations is primarily due to 
the decrease of the Hubble rate during inflation, which is 
common to both perturbations and yields n^ < 0. In the 
following analysis we also allow rii > 0, but we shall see 
that blue power spectra of this type are not favoured by 
the data. 

We then have four parameters, A — \J A^ + A?^ , -B, 
A, and n describing the effect of correlated perturba- 
tions, where n ~ \ ~\- 2ni is defined to coincide with the 
standard definition of the spectral index for adiabatic 
perturbations. We leave an investigation of the full six 
parameters for future work. 

By defining the entropy-to-adiabatic ratio B* = B/A 
the parameter A becomes an overall amplitude that can 
be marginalized analytically (see below). In the follow- 
ing, to simplify notation, we write A — I and drop the 
star from B* . We limit the analysis to B > and 
< A < TT, since there is complete symmetry under 
A ^ -A and under {B -^ ~B, A^tt-A). Further, we 
allow three background cosmological parameters to vary, 
LUb = ^bh"^ , Wc = ^cdmh"^ , and f^A where Vtb,cdm,A is 
the density parameter for baryons, CDM and the cosmo- 
logical constant, respectively. Since we assume spatial 
flatness, the Hubble constant is h? = """t,"'' . Our aim is 
therefore to constrain the six parameters 
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by comparison with CMB observations. We consider 
the COBE data analysed in |jl^, and the recent high- 
resolution Boomerang m, Maxima p] and DASI data 
B . In order to concentrate on the role of the primordial 
spectra (and limit the numerical computation required) 
we will fix the neutrino masses (zero) and spatial curva- 
ture (zero). We will also neglect any contribution from 
tensor (gravitational wave) perturbations. 

We use a CMBFAST code ^ modified in order to 
allow correlated perturbations to calculate the expected 
CMB angular power spectrum, Ci, for all parameter val- 
ues. (Our Ci is defined as d = l{l + l)C; /{2tt) where Cf 
is the square of the multipole amplitude) . The computa- 
tions required can be considerably reduced by expressing 
the spectrum for a generic value of B and A as a function 
of the spectra for other values. Let us denote the purely 
adiabatic and isocurvature spectra when i? = 1 as [C/]ad 
and [C;]iso respectively, and the correlation term for to- 
tally correlated perturbations i? = 1, A = as [C/]corr ■ 
Then we can write the generic spectrum for arbitrary B 
and A as 



Cl = [ClU + S2[C/]iso + 2BCOS A[G]e 
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We can obtain [C/]corr from Eq. (0) and using any 
ScosA 7^ 0. The library spectra [C;]ad and [C/]iso and 
[Cjjcorr can then be used to evaluate Ci for any B and 
A. A different set of library spectra will be needed for 



each set of cosmological parameters. When ni ^ 713 then 
A is not generally scale independent and so it would be 
necessary to evaluate the shape of the cross-correlation 
spectra [Cjjcorr for each form of A(fc), but one can always 
perform the scaling with respect to B analytically. 

The remaining input parameters requested by the 
CMBFAST code are set as follows: T^mh = 2.7267^, 
Yue = 0.24, N^ = 3.04. All our hkelihood functions be- 
low are obtained marginalizing over t^^ the optical depth 
to Thomson scattering, in the range (0,0.2) (larger t^ 
have a very small likelihood). We did not include the 
cross-correlation between band powers because it is not 
available, but it should be less than 10% according to [Q . 
An offset log-normal approximation to the band-power 
likelihood has been advocated by ||l^ and adopted by 
[Q,g|, but the quantities necessary for its evaluation are 
not available. Since the offset log-normal reduces to a 
log-normal in the limit of small noise we evaluated the 
log-normal likelihood 
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where Zg = log Ci, the subscripts t and d refer to the the- 
oretical quantity and to the real data, Cg are the spectra 
binned over some interval of multipoles centered on £i, ai 
are the experimental errors on Zg^d^ and the parameters 
are denoted collectively as aj . 

The overall amplitude parameter A can be integrated 
out analytically using a logarithmic measure dlogA in 
the likelihood. Analogously, we can marginalise over the 
relative calibration uncertainty of the Boomerang, Max- 
ima and DASI data (see [|l],p| ) , by an analytic integration 
to obtain the final likelihood function that we discuss in 
the following. We neglected beam and pointing errors, 
but we checked that the results do not change signifi- 
cantly even increasing the calibration errors by 50%. We 
assume a linear integration measure for all the other pa- 
rameters. 

In order to compare with the Boomerang, Maxima and 
DASI analyses we assume uniform priors as in |l|], with 
the parameters confined in the range i? e (0, 3), As 
(0,7r), n e (0.6,1.4), ujb G (0.0025,0.08), uj^ £ 
(0.05,0.4), r^A e (0,0.9). As extra priors, the value of h 
is confined in the range (0.45, 0.9) and the universe age 
is limited to > 10 Gyr as in |[. A grid of - 10,000 
multipole CMB spectra is used as a database over which 
we interpolate to produce the likelihood function. 

Figure 1 shows one of the best cases in our 
database, corresponding to (S, A,ri,a;fc, Wc^a) = 
(0.63, 7r/4, 0.9, 0.0225, 0.1, 0.7). The adiabatic ([G]ad), 
entropy {B^[Ci]iso) and correlated (2i3cos A[Ci]corr) 
components are shown. The primary effect of adding a 
positively correlated component is to reduce the height of 
the low-/ plateau relative to the acoustic peaks |jl^| . This 
is in contrast to the uncorrelated case where the addi- 
tion of entropy perturbations increases the plateau height 
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FIG. 1. Best-fit spectrum (solid line) and the component 
spectra, shown with the data with one sigma error bars, us- 
ing maximum-likelihood normalisation and unadjusted rela- 
tive calibration. 



relative to the peaks. Isocurvature perturbations only 
have a significant effect on intermediate angular scales 
for strongly blue-tilted spectra. They have a minimal ef- 
fect on the peak structure for n < 1. Thus we find a 
near-degeneracy between B and n when A = 0: the ef- 
fect of adding maximally correlated isocurvature pertur- 
bations mimics an increase in the primordial slope. This 
makes clear the importance of varying n when studying 
correlated isocurvature perturbations: a lower n allows a 
larger B to be consistent with the CMB data. 

This near-degeneracy is broken due to the effect of n 
on the slope at low-L In Fig. 2 we plot the likelihood for 
B and cos A, having marginalized over the other param- 
eters. The plot shows that the marginalized likelihood 
peak occurs for B = 0.4,cosA — 0.7, although the pure 
adiabatic case B = is well within one sigma. It is 
remarkable that when a non-zero correlation is allowed, 
quite large values of B become acceptable, up to _B = 1.5 
(to 95% c.l.) when cosA w 0.8. Anti-correlation, on 
the other hand, reduces the range of B. We also show 
the likelihood contours possible in a future Planck-like 
experiment with zero calibration uncertainty and accu- 
racy limited only by cosmic variance for I < 1000. This 
shows that future CMB data alone could detect a finite 
isocurvature contribution around the current peak of like- 
lihood. 

We found that the contour lines of the cosmological 
parameters wj, and lOc are almost parallel to B for B < 1. 
This means that the isocurvature perturbations do not 
alter significantly the best estimates for these cosmologi- 
cal parameters. On the other hand, increasing B enlarges 
the region of confidence for I^a and of n toward smaller 
values. 

Figure 3 summarizes our results: we plot the one- 
dimensional likelihood functions obtained by marginal- 
izing all the remaining parameters. Panel a shows that 
the contribution of isocurvature perturbations can be as 
large as the adiabatic perturbations, or even larger: we 
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FIG. 2. Contour plot of the two-dimensional likelihood for 
B and cosA. The contours enclose 40%, 86% and 99% of the 
likelihood and the star marks the peak. 

find that B < 1.3 to 95% c.l.. In contrast, uncorrelated 
isocurvature perturbations cannot exceed B < 0.5 to the 
same c.l.. The likelihood functions for n and J^a extend 
toward smaller values, as anticipated, while the CDM and 
the baryon density estimates remain quite unaffected. 
The average values are n = 0.94±0.1, ujb = 0.023±0.004, 
LUc = 0.1 ± 0.03, r^A = 0.72 ± 0.11, . 



unlike the case of independent perturbations where the 
relative height always decreases. Trotta et al ||l^ found 
that the CMB data was not consistent with a significant 
CDM isocurvature contribution because they restricted 
the primordial slope, n, to be unity. 

As can be seen from Fig. 3 our estimates of ojb and 
ujc are virtually unaffected by the addition of correlated 
CDM isocurvature perturbations. Thus, in our model, 
the nature of the isocurvature component can be inves- 
tigated almost independently of the composition of the 
matter component. 

The main conclusion of the present work is that the 
current CMB data is consistent with a large correlated 
CDM isocurvature perturbation contribution when the 
spectral slopes is allowed a tilt to the red (n < 1). The 
higher precision of future satellite data has the potential 
to detect the isocurvature contribution, if any, thereby 
showing that inflation was not a single-field process. 

The authors are grateful to David Langlois, Roy 
Maartens and Carlo Ungarelli for useful discussions. CG 
was supported by the ORS, and DW by the Royal Soci- 
ety. 
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FIG. 3. One-dimensional likelihood functions in arbitrary 
units. Green (light) dotted lines for the purely adiabatic mod- 
els (5 = 0); blue short-dashed lines for uncorrelated fluctua- 
tions (cos A = 0); red (dark) solid lines for correlated fluctu- 
ations. See text for further explanation. 

By contrast, Enqvist et al mm found that a large uncor- 
related isocurvature contribution is only consistent with 
blue tilted slopes. The reason for this difference is that 
correlations can cause the acoustic peak height to in- 
crease relative to the Sachs Wolfe plateau (see Fig. 1) 
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